Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN WF$t style file v2.2) 



The Luminosity Distance in Perturbed FLRW Spacetimes 



T. Pyne 

1101 East 58th Street, Chicago, IL 60615, USA 

M. Birkinshaw 

Department of Physics, University of Bristol, Tyndall Avenue, Bristol BS8 1TL 



Submitted: 9-July-2003 



ABSTRACT 

We derive an expression for the luminosity distance in FLRW spacetimes affected 
by scalar perturbations. Our expression is complete to linear order and is expressed 
entirely in terms of standard cosmological parameters and observational quantities. We 
illustrate the result by calculating the RMS scatter in the usual luminosity distance 
in flat (fimjty/O = (1-0,0.0) and (0.3,0.7) cosmologies. In both cases the scatter is 
appreciable at high redshifts, and rises above 11% at z = 2, where it may be the 
dominant noise term in the Hubble diagram based on SN la. 
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1 INTRODUCTION 

A substantial body of modern theoretical cosmology is 
concerned with the theory of small metric perturbations 
about Friedman-Lemaitre-Robertson- Walker (FLRW) cos- 
mologies. In particular, cosmologists often derive exact for- 
mulae for observables in such spacetimes and then apply 
standard statistical techniques in recognition of the stochas- 
tic nature of the perturbations. A great benefit of such a 
programme is that correlations between effects previously 
treated as distinct may be used to illuminate cosmological 
questions. 

As astronomers and cosmologists work with ever more 
distant sources, it becomes increasingly important to un- 
derstand and account for the deviations of our Universe 
from the simple FLRW spacetimes. One example, among 
many, where such an understanding is important is in the 
use of type la supernovae as standard candles to infer lumi- 
nosity distance as a function of source redshift, and hence 
to demonstrate that the Universe is currently in a phase 
of accelerated expansion (Reiss et al. 1998; Perlmutter et 
al. 1999). The inference that the Universe is accelerating 
was made on the basis of the luminosity distance formula 
appropriate for the background cosmology, not to the per- 
turbed cosmology in which we live. The gravitational effects 
of large-scale structure were taken to contribute to the er- 
ror budget at a level determined by the numerical studies 
of Wambsganss et al. (1997). Gravitational lensing by large- 
scale structure is known to produce effects on the order of 
those found in these numerical studies, and this is usually as- 
sumed to be the dominant effect on the luminosity distance 
(Kantowski, Vaughan & Branch 1995; Holz & Wald 1998; 



Holz 1998; Kaiser 1998). However, so far as we know, no 
rigorous exploration of this assumption has been attempted. 
In any case, the division of the effects of large-scale struc- 
ture into categories such as lensing, time-delays, and others 
is largely artificial: we measure luminosities and redshifts in 
the real Universe and no idealized FLRW background exists. 

In this paper we investigate the luminosity distance 
in linearly perturbed FLRW spacetimes. Our results are 
obtained by direct integration of the geodesic equations. 
A complementary approach has been employed by Sasaki 
(1987). Our results improve on those of Sasaki in two ways. 
First, we present explicit formulae for curved cases, whereas 
Sasaki's solutions in these cases require a prior integration 
of the perturbed geodesic equations. Second, our formulae 
are functions only of standard cosmological parameters and 
observational quantities. 

The plan of this paper is as follows. In section [5] we 
derive a formula for the luminosity distance, accurate to 
first order, in metric perturbed FLRW spacetimes. In sec- 
tion [3] we evaluate the cosmological weak lensing correction 
to the usual luminosity distance as a function of redshift in 
two model flat spacetimes, with fi m = 1 and (r2 m ,SlA) = 
(0.3,0.7). Section 0| summarizes our results. 

In this paper we use units such that G = c = 1. Greek 
indices /i, v, ... run over {0, 1, 2, 3}, and Roman indices i,j, ... 
run over {1,2,3}. The spacetime metric is taken to have 
signature +2, and the Riemann and Ricci tensor conventions 
are given by [V Q , V,a] u M = and R a g = R^a^p- 
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2 ANALYSIS 

In any spacetime there is a relation between the element 
of cross-sectional area of a radiation source, dA, and the 
solid angle that this source is observed to subtend at an 
observer, dQ. Although dA is unobservable, this relation is 
usually used to define a distance measure r, the angular- 
diameter distance, according to dA = r 2 dQ. In a general 
spacetime r will depend upon the particular locations of 
the source and observer, and possibly also on the choice of 
multiple null paths along which the source is observed. In 
perfectly homogeneous and isotropic FLRW spacetimes r 
depends only upon the observed redshift of the source and 
the cosmological parameters describing the evolution of the 
scale factor: the density parameters in matter, curvature, 
and cosmological constant and the Hubble constant, O m , 
Q K , Qa and Ho- In any spacetime, the luminosity distance, 
cZl, is related to r and the observed redshift z of the source 
according to dt = (1 + z) 2 r. In this section we exploit our 
previous solution (Pyne & Birkinshaw 1996) for r in linearly 
perturbed FLRW spacetimes to find dh- 

We choose the cosmological model to be an FLRW 
spacetime with embedded small scalar perturbations (f> (the 
extension to vector and tensor perturbations is conceptu- 
ally straightforward). The metric of such a spacetime may 
be written in the form 



ds 2 = a 2 [-(1 + 2<j))dr] 2 

+ (1 - 20) 7 -2 (dx 2 + dy + dz 2 )] 



(1) 



where 7 = 1 + ftf 2 /4, k is the spatial curvature parameter 
(±1 or 0), and f 2 = x 2 + y 2 + i 2 > 

Consider the radiation emitted by a source at event 
£, corresponding to conformal time rj e , that reaches an ob- 
server at event O at conformal time 770 at the spatial origin. 
Let the four-velocities of the source and observer be writ- 
ten < = (1/Oe) (l - </>e,«e) and < = (l/Oo) (l - <j>o, «o) 

respectively. 2 Suppose, further, that the observer chooses a 
spatial frame such that the source appears in the z-direction 
on the observer's sky. Then the relation between the cross- 
sectional area of the source at £ , dA, and the solid angle the 
source subtends at O, dQ,, was shown in Pyne & Birkinshaw 
(1996) to be 



dA = dQ (1 + 2Vo] 
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the integrals being taken along the path a;' ' = (77, 0, 0, r(ji)), 



1 A hat is used merely to distinguish this particular coordinate 
system, and not to denote a normalized vector. 

2 Note that these velocities are properly normalized if their spa- 
tial components are treated as first-order quantities. 



which is a null geodesic of the background spacetime, with 
affine parameter 77, provided that r = 2tan re ((77 — 77) /2). 
On this path 7 = sec 2 ((770 — rf)/2). In JHJl and the following, 
the subscript k on a trigonometric function denotes a set of 
three functions: for k = 1 the trigonometric function itself, 
for k — — 1 the corresponding hyperbolic function, and for 
k — the first term in the series expansion of the function. 3 
The terms in (|5J have the following interpretations: 

(special relativistic transformation \ 
of solid angle at O due to J 
observer's local peculiar velocity / 

^2„2 / metric factors in induced \ 

(1 — 2<f> c ) — I area two- form in 

\ rest frame of emitter / 

(effect of the gravitational field \ 
on the light rays as they travel I 
from source to observer / 



„,2 



In particular, the first two integrals in expression JHJ for 
det M describe longitudinal effects of the perturbations, 
while the last term describes the transverse effects usually 
interpreted as gravitational lensing. 

It is convenient to write in the form 
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dr\ (77 - 77c) <j> tV 



di] sin K (77 - 77e) (j> tV 



dr] cos K (77 - 77c) (j> 



dr] sin K (77 - 77c) 7=7 + <t>,yy] (4) 



which is achieved by noting that the directional derivative 
along x^°\ d(j>/dr] = (p^—'ytfii, and performing an integration 
by parts. 

While © is, in principle, a complete description of the 
angular-diameter distance in an FLRW spacetime perturbed 
by scalar perturbations, it suffers from a number of defects 
that limit its practical utility for observational cosmologists. 
Most seriously, it exhibits explicit dependence on the unob- 
servable point of emission £ rather than on the observable 
redshift. Less seriously, it does not explicitly depend on the 
usual cosmological parameters. We proceed to remedy both 
defects. 

Let Q at conformal time ?7 q and at spatial position 
(0, 0,f(?7q)) on a;' -' be an event with the property that the 
redshift of a source at Q observed at O and computed in 
the background metric is numerically the same as that of a 
source at £ observed at O computed in the full perturbed 
spacetime. Since 77c and 77 q will differ only at first order, and 
is valid only to this order, we can replace 77c by ?7 q in J1J 
without further modification. The metric two-form factors 
are trickier, however. Pyne & Birkinshaw (1996) showed how 
to calculate the coordinate displacements between £ and Q 



3 Pyne & Birkinshaw (1996) write, more precisely, A e where we 
write ?7e, but the two quantities differ only at first order, which 



renders them equivalent in (2) . 
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suffered by a photon due to the presence of the perturba- 
tions, (Sr), 5x l ). The components of the displacement needed 
here are those corresponding to conformal time and longi- 
tudinal coordinate z. They are 



Sz = 



"7q (ST) + I S ) 



(5) 
(6) 



where we have used a dot to denote a conformal-time deriva- 
tive. In the above, the perturbation to the time-like compo- 
nent of the photon wavevector 

n 

(7) 



fc (1)0 = -20 + 2 / d7/0,„ 
and 

rn<i 

Is = 2 dri4> 



is the usual expression for the Shapiro delay, both expres- 
sions being taken over the background path a;' '. 4 We point 
out that the displacement Sr] is easily interpreted as the 
sum of three effects: a differential doppler shift, a differen- 
tial gravitational redshift, and what is commonly termed the 
"integrated Sachs- Wolfe effect" in studies of the CMBR. 

The expressions 15161 allow us to Taylor expand the 
metric two-form factors about their values at Q with the 
result 5 
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Recognizing that along the background path r/7 = sin ft (?7o — 
77) we see 

^$ « ^ [l + 2^77 -2 cot K (7,0 -77^77 

7c 7q L flq 

-2cot K (77 -r? q )/s] • (10) 

Since (1 + z) 2 a q f q /^ q = d^\ the luminosity distance of the 
background cosmological model, equations (|5J, @ and 11 Oil 
lead to 



d L = d { °\l + S) 
with 



(11) 



8 = « + — Sr/ - cot K (r/ - T] q )Sri 

d q 

— cot K (f) — rj q )Is — 4>q — 2</>o 



In the notation of Pyne &: Birkinshaw (1996), 8rj 



SXa 



xWiXe) and Sx i = X^Ae) - X^^Aq) + X^^Xe). By Taylor 
expansion, 5x z fc q °^5A q + a;( 1 ' 1 (A c ). An integration by parts 
on the expression for zW^Ae) given in Pyne & Birkinshaw then 
produces 15161 . 

5 We note that the transverse derivatives in the Taylor expansion 
do not contribute once they are evaluated at Q. 
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dr) cos K (77 — r) q )(j) 



dr) sin K (77 — r/ q ) fj [4>^x + <t>,yy] (12) 



We now show how each term on the right-hand side of 
the above equation may be expressed in terms of observable 
quantities and standard cosmological parameters. First we 
change to the more familiar comoving coordinates (77, x, y, z) 
in which the background metric assumes the form 



/ , . dr 2 2 2 2 • ,, , 2 

+ r do + r sm 



a (77) drj 



(13) 



Here r = y^x 2 +y 2 + z 2 is the comoving radial coordinate 
distance from the origin. 6 This change is effected by the 
transformation x l = & /j. A number of useful results can 
now be derived. First, along x^°\ r — sin K (?7 — 77) and, 
in particular, sin K (77 — ?7 q ) = r q , the comoving coordinate 
distance to Q. From this follows cos K (?7o — 77) = yl — nr 2 
and the change in measure, dr — — cos K (r/ — rf)dr]. 

On x^°\ the Jacobian of the coordinate change may be 
shown to take the value 

h 




— = I i 



dxi 





1 1 ~1 A 



which secures, again on x 

<t>,xx = l 2 <t>xx 

4>,yy = l 2( t>yy 



(0) 



(14) 



(15) 
(16) 



Finally, it will be convenient to employ the common 
formulae for trigonometric functions of angular sums, in the 
forms 



sin K (77 - 77 q ) = sin K (77 - ?7 )cos K (770 - ?7 q ) 

+ cos K (77 - 77o) sin re (770 - 77 q ) 



cos K (77 - 77 q ) = cos K (77 - 770) cos K (770 - 77 q ) 

— k sin K (77 — r/ ) sin K (770 — 77 q ) 



= yl - Kr q yl - nr 2 — Krr q 

v-Vq. = -{vo - v) + (no - 



. -1 . . -1 

- sm K r + sm K r q 



(17) 

(18) 
(19) 



Using these identities, we find 
S = (v q — <f> — 2<f> q + hsw) 



dr 



-2r" 



vT 



sin K 1 r + sin K 1 r q ] <j} :V 



6 Some authors reserve this term for a r which is a proper dis- 
tance on the r] = r] hypersurface. r here is a dimensionless coor- 
dinate distance. 
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where 
-f/sw = — 2 



dr 



is the integrated Sachs- Wolfe effect, 
dr 

Is 



vl — Kr 2 
is the Shapiro effect and 



8r) = («q 



q + Iisw) 



(20) 
(21) 

(22) 
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Because of the way in which we defined Q, the ob- 
servable redshift of our source obeys the standard equation 
1 + z — a /a q and many of the standard cosmological equa- 
tions apply with their usual interpretation. In particular, 
the Hubble factor at the location of a source of redshift z is 
given by 7 



H = = H a E(z) 



with 



E{z) = + 2 ) 2 (1 + zfi m ) - 2(2 + z)Q A . 

Additionally, the comoving proper distance to Q is 

i«(vTO/i(*)) 



H oy /\n7\ 

where f\ (z) is the integral (shown in Fig. 
dl 



Mz) = 



o ^(0 



(24) 
(25) 
(26) 

(27) 



and we have the usual expressions for the density parameters 
in matter, vacuum energy, and curvature 



_ 8tt 



A 
3H* 



(28) 



respectively. As usual, these constants are related by the 
identity fi m + ^A+fi K = 1 which has the crucial consequence 
that for tt K / 



(29) 



and hence a is expressible in terms of the usual cosmological 
parameters. This is in contrast to the flat case where there 
is no meaningful curvature scale available. 

7 It is important to realize that the numerator here is a confor- 
mal time derivative. This accounts for the unusual power in the 
denominator. 





Figure 1. The functions fi{x) and f2(x) for (Q m ,nA) 
(1.0,0.0) (top) and (fi m ,n A ) = (0.3,0.7) (bottom). 



For the curved cases, then, both r q and a r q are func- 
tions of the standard cosmological parameters and the ob- 
served source redshift and so 11201 represents the desired 
formula for the luminosity distance. Simply for clarity we 
record, for these cases, 



ty, = sin. 



(30) 



and, valid for all cases, curved and flat, 



r q 8ri = i . '- (31) 



E{z) sin„ 



|n-l/i(*: 



In the fiat case our formula simplifies considerably, taking 
the form 

8 = (v* — 4> a — 20 q + Iisw) 
~r~ 1 Sr] - 2r~ 1 Is 



dr [— r + rj r [4> tXX + <f>, yy \ 



(32) 
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In contrast to the curved cases, in the flat case it is not 
possible to express r q in terms of standard cosmological pa- 
rameters and source redshift alone. As a result, it looks at 
first as if the terms in 132H involving integrations will prove 
problematic to our program. In fact, however, this difficulty 
is easily overcome by the change of variables r i— > r q w and 
k = |k| i— > v/r q performed after Fourier decomposition of <j> 
with mode-variable k. This procedure is illustrated in the 
examples below. 



3 LUMINOSITY DISTANCE SCATTER FROM 
WEAK LENSING 

The last term in 1321 . and its antecedent in 1201 . are pre- 
cisely the standard cosmological weak lensing terms usually 
analyzed in isolation. In this section we illustrate how these 
terms can be expressed as functions of the standard cosmo- 
logical parameters, and work out the contribution of the last 
term in l|32|l numerically for a simple f2 m = 1 flat cosmol- 
ogy, and for the currently popular flat (f2 m , £2 A ) = (0.3, 0.7) 
cosmology. The identification and analysis of the weak lens- 
ing term using a perturbed geodesic equation is not novel 
(Bernadeau, Van Waerbeke & Mellier 1997; Kaiser 1998). 
Its emergence as the dominant term in a physically sensible 
and complete description of the luminosity distance is, how- 
ever, and we include its analysis here for completeness as it 
is somewhat tricky to cast in the proper form. 
In order to evaluate 



Elfins — 



dr I — 



+ ■ 



(33) 



in either cosmological model it is necessary to specify the po- 
tential perturbation, <j>. <f> satisfies the perturbed Einstein's 
equations in the longitudinal gauge as given, for example, as 
equations (5.17)-(5.19) in Mukhanov et al. (1992) for f2 A = 
cosmologies. For the flat Q m — 1 cosmology during an epoch 
of matter-domination, and when the hydrodynamical per- 
turbations described by <j> are pressureless, <f> is independent 
of conformal time, taking the form 4>(x,r]) = <^>(x). 8 For 
models with vacuum energy, (f> varies with conformal time 
according to </>(x, rf) = s(ry)0(x) with s(rj) obtained from 



8{z) = {l + z)E{z)%& 
MO) 

where /2(a) is the integral 
(l + l)dl 



Hz) 



E{iy 



(34) 



(35) 



(sec Fig. and 2(77) is given through inversion of 

fx(z) = a Ho(r) - rf) . (36) 

The factor /2(0) in 134H ensures the conventional normaliza- 
tion for the Sl m = 1 case. 

While the temporal variation of the potential is deter- 
ministic in the models considered, the spatial variation is un- 
constrained. It is common to imagine the actual </>(x) as one 

8 It is possible to express the Laplacian of </>(x) in terms of the 
fractional perturbation to the hydrodynamical energy density on 
a surface of constant conformal time via what is usually termed 
a "cosmological Poisson equation" but this is not necessary here. 



realization drawn from a Gaussian random field. In this pic- 
ture, the random nature of <j> is usually expressed by specify- 
ing the probability distributions of its Fourier components, 



(x,r?) = s(r]) / h(k)e 



: j£k_ 

(2tt) 3 



with 

<h(k)fc*(k')> 
(ft(k)ft(k')) 
(h*(k)h*(k')) 



(27r) 3 <5(k - k')P(fc) 

(27r) 3 5(k + k')P(fc) 
(27r) 3 <5(k + k')P(fc) 



(37) 



(38) 
(39) 
(40) 



where the last two correlators follow from the first because 
of the requirement that <f> be real. 

Now we are ready to compute the RMS fractional devi- 
ation due to lensing. It is easiest to take the various corre- 
lators into account by writing cf> in the form 



(*,??) = 



-<Hx,j?) + ^(x,??)* 

s(n) 



= ^ / ft(k) 



' V 

s{rf) 
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: ri 3 k 
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whence 

(<5lens) 



(2 



J d 3 kP(fc) [kl + k. 



2 Y 



dr 



1 rs(r)e i 



(42) 



with s(r) given by the compositions of 13411 and 1361 above 
with the identity, true along ar ' where we need it, (r) ~rf) = 



3.1 Lensing effects from a simple power spectrum 

We must know the potential power spectrum P(k) to per- 
form the integral in 1421 . For the purposes of illustration, we 
follow Seljak (1994) and assume a broken power-law form 



P(k) = 



Ak~ 3 
Akik- 7 



for k < ko 
for k > ko 



(43) 



with ko, the turnover wavenumber, corresponding to a phys- 
ical scale of about 10 Mpc. The constant A is set by the 
CMBR quadrupole of Q\ = (6 X 10~ 6 ) 2 by 



Q 



2 2071-^1 



<) 



k 2 P(k)ji(2k/H )dk 



(44) 



(Bond & Efstathiou 1987). Here j2(x) is the spherical Bessel 
function of order 2 and K\ is the amplification coefficient of 
Kofman & Starobinskii (1985) which accounts for the time 
variation of the potential. For the f2 m = 1 case, K2 = 1 and 
for the (fi m ,fi A ) = (0.3,0.7) case K 2 = 1.19. The normal- 
ization for the power spectrum then becomes 



A 



11 

4.4 x lO" 11 



6.2 x 10" 11 if (f2 m ,n A ) = (1,0) 



if (n m ,n A ) = (0.3,0.7). 

We simplify the notation by generalizing 14311 to 
P(k) = Ak^- 3 k~ n 



(45) 



(46) 
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which can be used both above and below the turnover 
wavenumber, and change variables in the radial integration 
to w — r/r q and in the |k| integration to v 
this into 1421 . we obtain 



2.5 



(<5lens) — 



A(k r q 



(27T) 



dv v 



du(l 



dw w (w — 1) s(w, z) e 



(47) 



(48) 



for the contributions to (5^, ns ) at large and small k. In this 
equation, s(w,z) given by the composition of s(z), given in 
equation il'-> 1 1 . with z(w,z), obtained by solving 



wfi(z) = fi(S) 



(49) 



It may be worth noting here that z is simply an intermediate 
variable that will disappear from the final result whereas z 
is the observable redshift of the source. 

The source redshift and the cosmological parameters 
f2 m and Q.a can be seen to enter into (<5i 2 ena ) m the same 
manner: that is, via s(w, z) and the product vo = kor q which 
is equal to 2n times the comoving proper distance to Q in 
units of the turnover wavelength 



Vo 



kor„ = 



2tt 
Ao-ffo 



dl 

W) 



1884 ft" 1 fi(z) 



(50) 



If we use the values of n from 14311 . then the total lumi- 
nosity distance scatter {5f cns } is given by the explicit form 

hi 



(<5lens) 



A 
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dv v 



du (1 — u ) 



x / dww(w— 1) / dw'w'(w' — l) 



ii 



ii 



x s(w, z) s(w' , z) cos (uv(w — w')^ 
dv ' 

V 



A 4 



du (1 — u 2 ) 2 



1) 



< / dww(w— 1) / dw'w'(w'- 
o Jo 

Xs(w, z) s(w', z) cos (uv(w — it/)J . (51) 

While the calculation of the functions /i and /2 (Fig0, 
z (Fig. |2j, and s (Fig. present no difficulties, the four- 
dimensional integrals in I5H have some unpleasant charac- 
teristics, which can be brought out by further analysis. 

It is convenient to change the order of integration in 
151H . starting by performing the u- integral, which can be 
done analytically. Define 
l 

2\2 



g{x) = I dll 
48 



■ cos X + 



2 ) cos [ux) 
16 



(52) 



Then integral 151H can be rewritten 
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Figure 2. The function z(z,w) for (tt m ,n A ) = (0.3,0.7). Con- 
tours are drawn at z = 0.5, 1.0, . . ., 4.5. 




Figure 3. The function s(z,w) for (C m ,f2 A ) = (0.3,0.7). For 
(f2 m ,r2 A ) = (1.0,0.0), s = | for all redshift. Contours are drawn 
at s = 1.07, 1.10, ...,1.29. 



where the function 

el 



«(«) 



dyy g{yq) + 
ii ./ j 

i 



dy 



g{yq) 



dy [y 3 g (yq) + -g(- 
y \y 



(54) 



and q — vo (w — w'). t(q) is the key function in what follows, 
but once calculated, the same t(q) function can be used for 
all flat cosmologies, since the cosmological behaviour is con- 
cealed in the scaling quantity vo and in the (w,w r ) integra- 
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v -1 — 



series to q 
-a 

eries to q 



i°gio(q = To(w-W')) 

Figure 4. The function g 4 i(<?) for (S]„,f2 A ) = (0.3,0.7) and a 
range of q appropriate to the integral for redshift z = 2.5. Note the 
large range of t(q), the singularity as q — > 0, and the oscillations in 
t(q) at moderate and large q. This weighted q 4 t(q) has a maximum 
at q = 4.77: i(g) diverges as — lnq at small q. 



tions via the s{w,z) function. Fortunately, the integral in 
15411 can be performed, to yield 

16 



t{q) = — q 5 (— 36 sing + 6gcosg + 30g — 2q 2 sing 
to 



3 . 4 • 

j cos q + q sm q 



5 Ci(g)) 



(55) 



where the reader is to understand |g| for q (formally, the 
equation is slightly different for negative q, since t(q) is even 
in q but Ci(g) is defined with a cut along the negative q 
axis). The cosine integral function Ci(g) appearing in 15511 
causes the major difficulty in evaluating the integral 1531 1 
because as q — *• 0, 



*(9) 
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lng + -q 



' q 4 +0(q°)) (56) 



4032 



where 7 = 0.577... is Euler's constant. This exposes the 
logarithmic singularity at q — 0. By contrast, Ci(g) is well- 
behaved as g — » 00, with |Ci(g)| oc g _1 , so that in this limit 



t(q)~> 32 q~ 



1 sin g H 5- cos g + O 



(57) 



The form of t(q) is shown in Fig.|l] note that the asymp- 
totic forms I|56I570 are good descriptions of the behaviour of 
t(q), but that the logarithmic divergence at q — *• makes 
the integration in l|53jl improper, with implications for the 
numerical scheme adopted. 



3.2 Numerical results 

If we use 1431 for the potential power spectrum, then the 
calculation of the integral for (5f cns ) becomes a problem of 
covering the (w, w') plane adequately at w « w' so that 
the logarithmic singularity and the oscillations of t(q) for 
moderate w — w' are well sampled. An efficient solution is 
to alter the integration variables (w,w') to (w a ,iub), where 

w » = -i ( w ~ w ') 



(58) 



U'b 



(59) 



Since g = vow a , and is independent of Wi, the integration 
now needs to be done carefully only in the w a direction, while 
the change in the integrand in 153t in the Wh direction is 
slow. Furthermore, the symmetry (w,w') — ► (w',w) implies 
that only the triangle 
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(60) 



(61) 



needs to be included in the integral since the ™ a < triangle 
gives the same contribution. For the case (S7 m , Qa) = (1,0), 
a further simplification is possible, and the smaller triangle 

1 

71 



0, 



(62) 



(63) 



is all that is needed. 

The dashed line on Fig. shows the result obtained for 
(fimj^A) = (1,0). The effect rises approximately as 



v/<CJ=12xl0- 2 2 5 



(64) 



at small z (with h = 0.72), and reaches about 6% at z — 
1, and more than 10% at z = 2. The solid line on Fig. 
shows the result for (f2 m ,fiA) = (0.3,0.7). The curve has 
a similar shape to that for the Einstein-deSitter case, but 
larger lensing effects are seen at z > 1.7. At z < 0.4 this 
curve can be represented to an accuracy of 4% by 



Vkj = 7.5 x 10- 



(65) 



The results for both cosmological models show the z? de- 
pendence that one would have expected on the basis of 1511 
alone. Numerical results for both cases (accurate to better 
than 1%) are given in Tabled 



4 SUMMARY 

We have presented a complete formula 1201 for the luminos- 
ity distance in linearly perturbed FLRW spacetimes. The 
simpler form 132H is appropriate in flat spacetimes. These re- 
sults give the first explicit presentations of the various gravi- 
tational effects which modulate the unperturbed luminosity 
distance in terms that can be related to observable quanti- 
ties, as shown explicitly for (fi m ,fiA) = (1,0) and (0.3,0.7) 
cosmologies in Section f3. II 

The cosmological weak lensing term from I32H leads to 
a fractional scatter in the luminosity distance, (5i cns ), which 
can be significant as shown in Tableland Fig-El Clearly this 
effect is appreciable for quasars and high-redshift galaxies, 
especially as such objects are now being seen to 2 > 6. The 
principal assumptions required to derive these results are 
that the form I43H is a good description of the potential 
power spectrum, and that the time evolution is accurately 
described by 5(77). 

At z = 1, the value of y/{Sf cns ) is roughly 0.06. Ob- 
jects of fixed absolute magnitude would therefore show an 
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redshift 

Figure 5. The value of yj (<5[ 2 cns ), expressed as a percentage 
change, as a function of redshift to z = 2.5 for (Q m , Qa) = (1, 0) 
(dashed line) and (S7 m ,f2A) = (0.3,0.7) (solid line). 

apparent magnitude scatter of about 0.12 mag from lensing 
alone. Since the scatter of SN la absolute magnitudes (af- 
ter correction to a common light curve) is about 0.17 mag 
(Perlmutter et al. 1999), it is clear that lensing makes a 
significant, and increasingly important, contribution to the 
scatter of the SN la Hubble diagram as it is extended to 
redshifts > 1. At the redshift limit of the SNAP mission 
(z ~ 1.7; Perlmutter et al. 2003), the lensing- induced scat- 
ter of supernova apparent magnitudes rises about 0.2 mag, 
and becomes a dominant contribution to the intrinsic noise 
in the Hubble diagram. The associated Malmquist bias may 
also become important. 

While the usual cosmological weak lensing term is dom- 
inant, other terms are present and can be expected to af- 
fect the correlations between luminosity distance corrections 
and other physical quantities such as the integrated Sachs- 
Wolfe effect. The approach we have used, direct integration 
of the null geodesic equation, can also be used to examine 
the effects of vector and tensor perturbations, and can be 
extended, with somewhat more difficulty, to higher orders 
using the results of Pyne & Carroll (1996). 



Table 1. Results for yj (<5[ 2 ens ) for the two cosmologies 
Redshift, z ^(5^) (%) for (fi m , Q A ) = 





(1.0,0.0) 


(0.3,0.7) 


0.01 


0.014 


0.007 


0.02 


0.040 


0.021 


0.03 


0.073 


0.038 


0.04 


0.11 


0.059 


0.05 


0.15 


0.082 


0.06 


0.20 


0.11 


0.07 


0.25 


0.14 


0.08 


0.30 


0.17 


0.09 


0.35 


0.20 


0.10 


0.41 


0.23 


0.20 


1.05 


0.64 


0.30 


1.77 


1.16 


0.40 


2.50 


1.74 


0.50 


3.22 


2.37 


0.60 


3.93 


3.02 


0.70 


4.61 


3.70 


0.80 


5.26 


4.37 


0.90 


5.89 


5.05 


1.00 


6.49 


5.73 


1.20 


7.62 


7.04 


1.40 


8.64 


8.31 


1.60 


9.59 


9.51 


1.80 


10.5 


10.6 


2.00 


11.3 


11.7 


2.50 


13.0 


14.0 


3.00 


14.5 


16.2 


3.50 


15.7 


18.0 


4.00 


16.8 


19.6 


4.50 


17.8 


21.0 


5.00 


18.6 


22.3 
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